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Abstract 

The recent results of Forshaw and Sabio Vera on small- x jet rates to order of are 
extended to all orders, for any number of jets. A simple generating function is obtained. 
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1 Introduction 



The summation of logarithms of 1/x in deep inelastic structure functions at small values of 
Bjorken x leads to the Balitskii-Fadin-Kuraev-Lipatov (BFKL) equation |], [|, which in the 
leading approximation sums terms of order [a s ln(l/x)] n . Recently the next-to-leading terms 
have also been computed H [| . 

The usual derivation of the BFKL equation assumes the dominance of multi-Regge kine- 
matics (i.e. strong ordering in the Sudakov variables). This is valid for the calculation of the 
totally inclusive structure functions, in which one sums over all hadronic final states. When 
studying the structure of the final states, however, one should take account of QCD coherence, 
which effectively imposes an angular ordering constraint on the emission of soft gluons || |7]]. 
The resulting 'CCFM' formulation [H reduces to the BFKL equation for the inclusive struc- 
ture functions, but leads in general to different exclusive multi-gluon distributions from those 
expected from multi-Regge kinematics || . 

In a recent paper, Forshaw and Sabio Vera showed that, in the leading logx (LLx) 
approximation, to third order in a s , the rates for emission of fixed numbers of 'resolved' final- 
state gluons, together with any number of unresolvable ones, are the same in the multi-Regge 
(BFKL) and coherent (CCFM) approaches. Here 'resolved' means having a transverse momen- 
tum larger than some fixed value /i R , and the LLx approximation means keeping only terms 
that have two large logarithms for each power of a s , at least one of which is ln(l/x). In this 
approximation, each resolved gluon can be equated to a single jet, since to resolve it into more 
than one jet would cost extra powers of a s with no corresponding powers of ln(l/x). 

The present paper extends the work of Forshaw and Sabio Vera to all orders, for any number 
of resolved gluons. The BFKL and CCFM formulations are shown to give the same jet rates in 
LLa; approximation to all orders. The factorization of collinear singularities is demonstrated, 
and a simple generating function for the jet multiplicity distribution is obtained. 



2 Multi-Regge (BFKL) analysis 



As in Ref. 0, we start from the unintegrated structure function of a single gluon, f(x, k), which 
in the exclusive form of the BFKL approach satisfies the equation 

f(x, k) = 6(1 - x) 5 2 (k) + a s J ^l^A R (z, k)6(q - fi) f(x/z, q + k). (1) 

Here a s = 3a s /n, k is the transverse momentum of the gluon probed in the deep inelastic 
scattering, q is that of an emitted gluon, fi is a collinear cutoff (which cancels in the inclusive 
treatment) and Ar is the Regge form factor 

( I k 2 \ 

A R (z, k) = exp ( -a s ln-ln— j . (2) 



The contribution from emission of n gluons is obtained by iteration, 

/<">(*, k) = f[f ^^a s A R (z h k t )9( gi - p)6(x - x n )5\k /,, ) . 



(3) 



where 



Xi = J] Zj , fci = - ^ ?j . (4) 

The contribution to the structure function at scale Q is then obtained by integrating over all 
H <qi <Q: 



FV>(x,Q,ri = H[ 



dqf d(f>i dzi 



i J f l Qi 



qf 2tt Zi 



a s A R (zi,ki)S(x - x n ) . 



(5) 



2.1 One-jet rate 

Now consider the effect of requiring one emitted gluon, say the jth, to have qj > /i R while 
all the others have < /i R . This defines the contribution of one resolved gluon (plus n — 1 
unresolved), F( n ' ljet ): 



" [»l dqf dfa dzi _ 

Notice that for i < j the contribution is identical to the (j — l)-gluon contribution to the 
structure function evaluated at x — Xj-i and Q = /i R . On the other hand for i > j we have 
ki = qj in leading logarithmic approximation, and so the qi integrations become trivial: 



F^\x,Q,^,im) = Ysfdx^F^ix^^tif 

• =1 Jo J fi 



Q 2 dqUzj . , 

-±—a s A R (zj, qj ) 

Mr Qj z j 



n 



dzi 



2a s SA R (zi, qj)5(x - x n ) 



(7) 



i=j+i • 

where S = ln(/x R //x). 

To carry out the integrations it is convenient to use a Mellin representation, 

r-l 



F„ 



dxx UJ F(x) , 



(8) 



so that 



Zb j = l J ^R Qj 

rQ 2 dq) 
V| ~q] 



j 



—^A R (z,qj) 

z 



z 

2a^S 



n-j+l 



-1 n-j+l 



uj + a s ln(q]/fj, 2 ) 



(9) 



FM*\Q, /i R , fi) = F w (// a> fi) l Q ^H^q, /i R ) 



Summing over all j and n gives the total one-jet contribution. 

fQ dq 
'm q 

where 
and hence 



2a s 

a; + 2a s ln(g/ /i F 



FpMiQ,^,!*) = F^n) In (l + -^t) 
where T = ln(Q//i R ). The asymptotic behaviour of the structure function is 



F w (// R ,/i)~ f-|J = exp[2S , 7(a s ,cu)] 



where 7(a s ,u;) is the Lipatov anomalous dimension: 



7 (a s ,^) = ^ + 2C(3)P) + ... 



(10) 

(11) 
(12) 

(13) 
(14) 



Taking the leading term in 7 gives the result of Ref . || , extended to all orders: 
i^ ljCt) (Q, Ma, A 4 ) = exp (^5) In (l + 



2«c 



T 



2a c 



TS - -T 2 
2 



2a c 



UJ 



-T 3 - -T 2 S + -TS 2 
3 2 2 



+ 



/2as 



UJ 



-TS 3 - -T 2 S 2 + -T 3 S - -T 4 
6 4 3 4 



+ 



(15) 



Notice that the collinear-divergent part (the S- dependence) factorizes out, and the fraction of 
events with one jet is given by the cutoff-independent function 

R^ et \Q,^) = ^pfq^^ = ^ (l + ^T) exp[-2T 7 (a s ,c)] . (16) 
2.2 Multi-jet rates 

Now suppose we resolve r gluons with transverse momenta qj > fi R . To leading logarithmic 
accuracy, the Regge form factors beyond the first of these have hi fixed at the largest of the 
q/s resolved so far, and therefore Eq. flip]) becomes 



r rQ da- 

j=i J w. Qj 



(17) 



where kj = maxj<j{gi}. Introducing t = ln(g//i R ), we have the general problem of evaluating 

r r T 

(18) 



G { :\T) = II / H u (m a x{U}) dtj . 

ZJy JO V / 



Defining Gffl (T) = 1 and introducing the generating function 

oo 

GU«,T) = X> r GjW(r), 

we have 



(19) 



r=0 



Gw(u,T) = exp 
= exp 

From Eq. (|TT|) we find in this case 

2a s 



dt (uH u (t) + uHHlit) + uH 2 Hl(t) + 



o l-utH w (t) 



(20) 



+ 2a s t 

Thus the r-jet rate is given by 



G u (u,T) 



i + (i-„)^r 



(21) 



where the jet-rate generating function R w is given by 



(22) 



u=0 



Ru{u,T) = exp f -— 



1 + (1 - u) T 

uj 



(23) 



The remarkably simple expression fl2~3|) is the main result of the present paper. We shall 
see that the same result is obtained from the CCFM formulation of small-x dynamics. After 
convolution with the measured gluon structure function, it gives the predicted jet rates in the 
LLx region ln(l/x) ^> T = ln(Q/yU R ) 1, to all orders in a s . 

Expanding to fourth order, we have explicitly 
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(24) 



Note that these sum to unity as expected. 
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From the generating function ( p3|) we can deduce other interesting quantities to all orders, 
for example the mean number of jets, 



<r>= JU,(u,!T) 



2a* 1 f2a s \ 2 . , 

T + - — T , (25) 



, u 2 \ w 



and the mean square fluctuation in this number, 

< r >"« =17 r+ 2(-7 T ) + 3(- r ) • < 26 > 

In general, the pth central moment of the jet multiplicity distribution is a polynomial in a s T/uj 
of degree 2p — 1, indicating that the distribution becomes relatively narrow in the limit of very 
small x and large Q/// R . 



3 Analysis including coherence (CCFM) 

Taking account of QCD coherence gives the angular ordering constraint qi > z^iq^i and in 
place of Eq. (||) one obtains 

" fQ 2 dq 2 dd>i dzi 

F { - n \x,Q,ix) = 11/ - « S A(^, fcj,fc)0(ft - Zi_igi_i)5(x - x n ) , (27) 
i=1 Jo gf 27T z s 

where we set z q = fi and A is the 'non-Sudakov' form factor 

( 1 k 2 \ 

A(z,k,q) =exp f -a s ln-ln — J . (28) 

Corresponding to Eq. (0) we now have for one resolved emission 

F^ et )(*,Q,/i R ,/i) = F^~ 1 \xj-\, /i R , /i) £ /^r^^.^.^ 

"'Mr 5,' JO 2j 



- r^R dq- dzi 

II / ~ T — -«sA(zi, g,-, g»)5(a; - x n ) . (29) 

Performing the Mellin transformation and summing over j and n, we obtain in place of Eq. fllQ|) 

Fi 1 J et )(Q,/i R , / i) = F L ,(/i R ,/i) / Q ^(g,g,/i R ) (30) 

g 

where 

r^n/l dz 

H w (k,q,/i n ) = 2a s — z UJ A(z,k,q)K UJ (z,k,q,fi B ), (31) 

JO -2 

Ka,(z, k,q,fi n ) = l + 2a s [ -^-^ r z' UJ A(z',k,q')K UJ (z\k,q',n R ). (32) 

Jzg g £ 
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The CCFM treatment predicts the same behaviour as BFKL for the structure function, and 
so the factor F u (p K , /i) in Eq. (|30D still has the form (|T3|). The explicit solution of Eq. (|3^ ) is 

K^z, fc, g, /i R = — — - . 33 
w + 2a s ln(fc//i R ) 

Introducing s = ln(fc//i H ), t = ln(g/// R ), we can write the expression in Eq. fl31| ) as 

H^s,t)= 2as _ exp (a a (t 2 - 2st) - cot) . (34) 

The exponential factor contributes only sub-leading corrections, i.e. terms with fewer factors 
of ln(l/a;) than of a s , which we are neglecting. Substituting in Eq. (|30|), we therefore find 
that to the same precision the CCFM result for the one-jet contribution is equal to the BFKL 
prediction (TO). 

To find the multi-jet fractions we use Eq. (|T9|) with 



G { :\T) = f[ f T H w fmax{^}, t 3 ) dt 3 , (35) 
j=l ^° ^ — J ' 

i^c now being a function of two variables, given by Eq. (0). However, H w depends only on 
its first argument, apart from the negligible exponential factor. Therefore the CCFM rates for 
more than one jet also have the BFKL values, given by Eq. (|23|). This completes the all-orders 
extension of the results of Ref . 101 . 



4 Conclusions 

The above results show that the multi-jet rates in deep inelastic scattering at small x, as defined 
here, are sufficiently inclusive quantities to be insensitive to the differences between the BFKL 
and CCFM formulations of small-x dynamics at the leading-logarithmic level. In both cases 
they are given by the simple generating function (|23|) . Differences would be expected at the 



sub-leading level and in more differential quantities such as multi-jet rapidity correlations |10 

Bearing in mind the importance of sub-leading corrections to structure functions at small x 
[0, [|, one would not expect these result to be of direct phenomenological relevance, although 
they may provide a useful cross-check on the results of numerical simulations of small-x final 
states 0, HI, [L3|. 
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